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Geometric Adaptive Control for Aerial Transportation of a Rigid Body 

Taeyoung Lee 


Abstract —This paper is focused on tracking control for a 
rigid body payload, that is connected to an arbitrary number 
of quadrotor unmanned aerial vehicles via rigid links. A 
geometric adaptive controller is constructed such that the 
payload asymptotically follows a given desired trajectory for 
its position and attitude in the presence of uncertainties. 
The coupled dynamics between the rigid body payload, links, 
and quadrotors are explicitly incorporated into control system 
design and stability analysis. These are developed directly on the 
nonlinear configuration manifold in a coordinate-free fashion 
to avoid singularities and complexities that are associated with 
local parameterizations. 

I. Introduction 

By utilizing the high thrust-to-weight ratio, quadrotor 
unmanned aerial vehicles have been envisaged for aerial 
load transportation [1], [2], [3]. Most of the existing results 
for the control of quadrotors to transport a cable-suspended 
payload are based on the assumption that the dynamics of 
the payload is decoupled from the dynamics of quadrotors. 
For example, the effects of the payload are considered as 
arbitrary external force and torque exerted to quadrotors [2]. 
As such, these results may not be suitable for agile load 
transportation where the motion of cable and payload should 
be actively suppressed. 

Recently, the full dynamic model for an arbitrary number 
of quadrotors transporting a payload are developed, and 
based on that, geometric tracking controllers are constructed 
in an intrinsic fashion. In particular, autonomous transporta¬ 
tion of a point mass connected to quadrotors via rigid links 
is developed in [4]. It has been generalized into a more 
realistic dynamic model that considers the deformation of 
cables in [5], and also the attitude dynamics of a payload, 
that is considered as a rigid body instead of a point mass, 
is incorporated in [6]. However, these results are based on 
the assumption that the exact properties of the quadrotors 
and the payload are available, and that there are no external 
disturbances, thereby making it challenging to implement 
those results in actual hardware systems. 

The objective of this paper is to construct a control system 
for an arbitrary number of quadrotors connected to a rigid 
body payload via rigid links with explicit consideration 
on uncertainties. A coordinate-free form of the equations 
of motion that have been developed in [6] is extended 
to include the effects of unknown, but fixed forces and 
moments acting on each of the quadrotors, the cables, and 
the payload. A geometric nonlinear adaptive control system 
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is designed such that both the position and the attitude of 
the payload asymptotically follow their desired trajectories, 
while maintaining a certain formation of quadrotors relative 
to the payload. 

The unique property is that the coupled dynamics of the 
payload, the cables, and quadrotors are explicitly incorpo¬ 
rated in control system design for agile load transportations 
where the motion of the payload relative to the quadrotors 
are excited nontrivially. Another distinct feature is that the 
equations of motion and the control systems are developed 
directly on the nonlinear configuration manifold intrinsically. 
Therefore, singularities of local parameterization are com¬ 
pletely avoided. 

As such, the proposed control system is particularly useful 
for rapid and safe payload transportation in complex terrain, 
where the position and attitude of the payload should be 
controlled concurrently. Most of the existing control systems 
of aerial load transportation suffer from limited agility as 
they are based on reactive assumptions that ignore the inher¬ 
ent complexities in the dynamics of aerial load transporta¬ 
tion. The proposed control system explicitly integrates the 
comprehensive dynamic characteristics to achieve extreme 
maneuverability in aerial load transportation. To the author’s 
best knowledge, nonlinear adaptive tracking controls of a 
cable-suspended rigid body with uncertainties have not been 
studied as mathematically rigorously as presented in this 
paper. 

II. Problem Formulation 

Consider n quadrotor UAVs that are connected to a 
payload, that is modeled as a rigid body, via massless links 
(see Figure [T]). Throughout this paper, the variables related to 
the payload is denoted by the subscript 0, and the variables 
for the i- th quadrotor are denoted by the subscript i, which is 
assumed to be an element of 1 = {1, • • • n} if not specified. 
We choose an inertial reference frame {ei, e* 2 , es} and body- 
fixed frames {bj 1 , bj 2 ,bj 3 } for 0 < j < n as follows. For the 
inertial frame, the third axis e 3 points downward along the 
gravity and the other axes are chosen to form an orthonormal 
frame. 

The location of the mass center of the payload is denoted 
by Xo G M 3 , and its attitude is given by R 0 G SO(3), where 
the special orthogonal group is defined by SO(3) = {R G 
M 3x3 | R t R = /, det[ii] = 1}. Let pi G M 3 be the point 
on the payload where the i-th link is attached, and it is 
represented with respect to the zeroth body-fixed frame. The 
other end of the link is attached to the mass center of the 
i-th quadrotor. The direction of the link from the mass center 
of the i-th quadrotor toward the payload is defined by the 
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Fig. 1. Dynamics model: n quadrotors are connect to a rigid body mo 
via massless links li. The configuration manifold is M 3 x SO(3) x (S 2 x 
SO(3)) n . 

unit-vector qi G S 2 , where S 2 = {q G M 3 | \\q\\ = 1}, and 
the length of the i-th link is denoted by h G M. 

Let xi G M 3 be the location of the mass center of the 2 -th 
quadrotor with respect to the inertial frame. As the link is 
assumed to be rigid, we have Xi = x o + Ropi — hqi- The at¬ 
titude of the 2 -th quadrotor is defined by Ri G SO(3), which 
represents the linear transformation of the representation of 
a vector from the 2 -th body-fixed frame to the inertial frame. 

In summary, the configuration of the presented system 
is described by the position xo and the attitude Ro of the 
payload, the direction qi of the links, and the attitudes Ri of 
the quadrotors. The corresponding configuration manifold of 
this system is Q = M 3 x SO(3) x (S 2 x SO(3)) n . 

The mass and the inertia matrix of the payload are denoted 
by mo G R and Jo G M 3x3 , respectively. The dynamic model 
of each quadrotor is identical to [7]. The mass and the inertia 
matrix of the 2 -th quadrotor are denoted by rrii G M and fa G 
M 3x3 , respectively. The i -th quadrotor can generates a thrust 
—fiRie 3 G M 3 with respect to the inertial frame, where fa G 
M is the total thrust magnitude and e 3 = [0, 0, 1] T G M 3 . It 
also generates a moment Mi G M 3 with respect to its body- 
fixed frame. The control input of this system corresponds to 
{fa, Mi}i<i< n - 

In this paper, the external disturbances are modeled as 
follows. The disturbance force and moment acting on the 
payload, namely A^, A Rq G M 3 are expressed as 

A* 0 = $ Xo (t,q,c\)0 Xo , A Rq = $ Ro (t,q,q)0 Ro , (1) 

where ^> Xo ,T>^ 0 : R x TQ R 3xn e denote matrix¬ 
valued, known function of the time t and the tangent vector 
(q,q) G T q Q of the configuration manifold, i.e., q = 
(x 0 ,R 0 ,qi,...,q n ,Ri,...R n ), and 0 Xo ,0 Ro e R n(,xl are 
fixed, unknown parameters for some uq. This type of uncer¬ 
tainties are popular in the literature of adaptive controls, and 
they may represent various modeling errors or disturbances, 
such as the uncertainties in the mass and the inertia matrix 
of the payload. Similarly, the disturbance force and moment 
acting on the 2 -th quadrotors are given by 

A*. = <& x . (f, q, q)0 Xi , A r . = & Ri (t,q,q)0 Ri , (2) 


where & x .,<f> R . : R x TQ —M 3xn0 and 0 Xi ,0 Ri G 
R n0 xi. Here, the disturbance forces are represented with 
respect to the inertial frame, and the disturbance moments 
are represented with respect to the corresponding body-fixed 
frame. 

Throughout this paper, the 2-norm of a matrix A is 
denoted by ||A||, and its maximum eigenvalue and minimum 
eigenvalues are denoted by A m [A\ and A m [A ], respectively. 
The standard dot product is denoted by x • y = x T y for any 

x, y G M 3 . 

A. Equations of Motion 

The kinematic equations for the payload, quadrotors, and 
links are given by 

qi=UiX qi = Qjiq u ( 3 ) 

Ro = Ro&o, Ri = Ri&i, ( 4 ) 

where uji G M 3 is the angular velocity of the 2 -th link, 
satisfying qi • ui = 0, and fiq and &i G M 3 are the angular 
velocities of the payload and the 2 -th quadrotor expressed 
with respect to its body-fixed frame, respectively. The hat 
map • : M 3 so (3) is defined by the condition that 
xy = x x y for all x, y G M 3 , and the inverse of the hat 
map is denoted by the vee map V : so(3) M 3 , where 

so(3) denotes the set of 3 x 3 skew-symmetric matrices, i.e., 
so(3) = {S G M 3x3 | S T = —S}, and it corresponds to the 
Lie algebra of SO(3). 

We derive equations of motion according to Lagrangian 
mechanics. The velocity of the 2 -th quadrotor is given by 
Xi = x o + Ropi ~ hqi- The kinetic energy of the system 
is composed of the translational kinetic energy and the 
rotational kinetic energy of the payload and quadrotors: 

T = ^TTT-o ||*o II 2 + ^0 • JqSIq 

n ^ 2 . 1 

+ ^ ^ -777^11 Jo H" R()Pi T" 2^* ' 

i= 1 

The gravitational potential energy is given by 

n 

U = -m 0 ge 3 ■xq-'^ m id e 3 • (^o + Ropi ~ kQi), (6) 

i =1 

where the unit-vector e 3 points downward along the grav¬ 
itational acceleration as shown at Fig. [T] The resulting 
Lagrangian of the system is C = T — U. 

The corresponding Euler-Lagrange equations have been 
developed according to Hamilton’s principle in [6], Here, 
it is generalized to include the effects of disturbances via 
the Lagrange-d’Alembert principle. Let the action integral 
be 0 = f to f C dt. Next, let the total control thrust at the 
2 -th quadrotor with respect to the inertial frame be denoted 
by Ui = —faRiCs G M 3 and the total control moment at 
the 2 -th quadrotor is defined as Mi G M 3 . There exist the 
disturbances A Xo , A Rq for the payload, and the disturbances 
A Xi , A Ri for the 2 -th quadrotor. The virtual work can be 




written as 


rt / n 

5W= ' {^o + RofjoPi - x qi} 

Jt 0 i=1 

n 

+Ew + A R.) • r]i + A Xo • Sx o + A Ro • r/o dt. 

i= 1 

The Lagrange-d’Alembert principle states that $0 = — 5W 
for any variation of trajectories with fixed end points. This 
yields the following equations of motion (see [ 8 ] for detailed 
derivations), 

n 

M q (x 0 - ge 3 ) -’Y^m i q i qlRopitlo = A Xo 

i= 1 
n 

+ + A Ii - m ih\\ui\\ 2 qi - rniqiqf Ro&lpi, (7) 

7=1 

n 

(J 0 - ^ miPiR^qjcftR 0 pi)fl 0 

7=1 

n 

+ ^ rriiPiRo q t qj (x 0 - ge 3 ) + O 0 Jo^o = A flo 

i=l 

n 

+ y^/5ii?o ( M i + A|. - - mtqiqjR 0 (llpi), 

i =1 

( 8 ) 

A; = -j-qi(x 0 - 5e 3 - RopMo + iiofigpi) 

H 

--7- ®(“i* + A *,)> ( 9 ) 

Ji&i + x = Mi + A#., (10) 

where = m y I + ^ M 3x3 , which is 

symmetric, positive-definite for any qi. 

Recall the vector xx^ G M 3 represents the control force at 
the i-th quadrotor, i.e., Ui = —/^e 3 . The vectors 77 - and 
uf- G M 3 denote the orthogonal projection of xx^ along qi , 
and the orthogonal projection of xx^ to the plane normal to 
qi , respectively, i.e., 

“I 1 = qiqfu u (li) 

ut = -qjut = (/ - qiqj)ui. ( 12 ) 


Therefore, u t = v)j + up Throughout this paper, the 
subscripts || and 1 of a vector denote the component of the 
vector that is parallel to qi and the other component of the 
vector that is perpendicular to Similarly, the disturbance 
force at the 7 -th quadrotor is decomposed as 

A l = q.q^xA, = ^i0 Xi , (13) 

a i = {i- qi qT)* Xi e Xi = *i8 Xi - (14) 


B. Tracking Problem 

Define a fixed matrix ? G M 6x3n as 


V = 


hx3 

Pi 


hx3 

Pn 


(15) 


Recall that pi describe the point on the payload where the 
x-th link is attached. Assume the links are attached to the 
payload such that 


rankp] > 6 . 


(16) 


This is to guarantee that there exist enough degrees of 
freedom in control inputs for both the translational motion 
and the rotational maneuver of the payload. The assumption 
( T 6 | ) requires that the number of quadrotor is at least three, 
i.e., n > 3. 


It is also assumed that the bounds of the disturbance forces 
and moments are available, i.e., for known positive constant 
B<$> , Be G M, we have 


(17) 

maxilla ||,||6» flo ||,|| [6 xt ||,..., \\9 Xn ||, 

(18) 


Suppose that the desired trajectories for the position and 
the attitude of the payload are given as smooth functions of 
time, namely xo d (t) G M 3 and Ro d (t) G SO(3). From the 
attitude kinematics equation, we have 

Ro d (t) = Ro d (t)(l 0d (t), 

where flo d (t) £ R 3 corresponds to the desired angular 
velocity of the payload. It is assumed that the velocity and 
the acceleration of the desired trajectories are bounded by 
known constants. 

We wish to design a control input of each quadrotor 
{/i,Mi}i<K n such that the tracking errors asymptotically 
converge to zero along the solution of the controlled dynam¬ 
ics. 


III. Control System Design For Simplified 
Dynamic Model 

In this section, we consider a simplified dynamic model 
where the attitude dynamics of each quadrotor is ignored, 
and we design a control input by assuming that the thrust 
at each quadrotor, namely Ui can be arbitrarily chosen. It 
corresponds to the case where every quadrotor is replaced 
by a fully actuated aerial vehicle that can generates a thrust 
along any direction arbitrarily. The effects of the attitude 
dynamics of quadrotors will be incorporated in the next 
section. 

In the simplified dynamic model given by 0 - 0 . the 
dynamics of the payload are affected by the parallel com¬ 
ponents u\ of the thrusts, and the dynamics of the links are 
directly affected by the normal components u t of the thrusts. 
This structure motivates the following control system design 
procedure: first, the parallel components u\ are chosen such 
that the payload follows the desired position and attitude 
trajectory while yielding the desired direction of each link, 
namely qi d G S 2 ; next, the normal components uf- are 
designed such that the actual direction of the links qi follows 
the desired direction q id . 






A. Design of Parallel Components 

Let di G M 3 be the acceleration of the point on the payload 
where the i -th link is attached, that is measured relative to 
the gravitational acceleration: 

di= x 0 - ge 3 + Ro&lpi - Rafail 0 . (19) 

The parallel component of the control input is chosen as 

u\ = pi + mili\\uji\\ 2 qi + rriiqiqf ai , (20) 

where pi E M 3 is a virtual control input that is designed 
later, with a constraint that pi is parallel to <^. Note that the 
expression of is guaranteed to be parallel to q L due to the 
projection operator qiqj at the last term of the right-hand 
side of the above expression. 


valid in general, as each pi is constrained to be parallel to 
qi. Instead, we choose the desired value of pi , without any 
constraint, such that 


^2lM d = F d , ^2 PiRo Pi d = M d , 


(25) 




i= 1 


or equivalently, using the matrix V defined at © 


~RoPu 




Rq Fd 

_Rq Pn d _ 


. M d . 


The motivation for the proposed parallel components be- 
comes clear if (20) is substituted into ([7|-([8]) and rearranged 
to obtain 

~Pu 

= diag[i? 0 , • • • i? 0 ] V T (VV T )- 1 

R-0 Rd 

. M d . 

n 

_Pn d _ 




m 0 (x 0 — ge 3) — A Xo + ^^(Pi + A|.), (21) 


i= 1 


Jotio + n 0 J 0 n 0 — Ar q + ^2 PiRo (Pi + A|.). ( 22 ) 

Therefore, considering a free-body diagram of the payload, 
the virtual control input pi corresponds to the force exerted 
to the payload by the i- link, or the tension of the i-th link 
in the absence of disturbances. 

Next, we determine the virtual control input pi. As in [9], 
define position, attitude, and angular velocity tracking error 
vectors e XQ , e# 0 , e^ 0 G M 3 for the payload as 


From the assumption stated at ( [T6] ), there exists at least one 
solution to the above matrix equation for any Fd, Md. Here, 
we find the minimum-norm solution given by 


(26) 


The virtual control input pi is selected as the projection of 
its desired value pi d along q i9 

Pi = (Pi d ■ Qi)Qi = QiQi Pi d i ( 27 ) 

and the desired direction of each link, namely qi d G S 2 is 
defined as 


Qu = 


Pi d 
II P'id I 


(28) 


It is straightforward to verify that when qi — q id , the resultant 
force and moment acting on the payload become identical to 
their desired values. 


&xo — Xo %0d 5 

e R 0 = 2 - 

£n 0 = — RoRo d n 0d . 

The desired resultant control force Fd G M 3 and moment 
Md G M 3 acting on the payload are given as 


Rd — m 0 ( k XQ e XQ ^xq^xq H - %o d ge 3) 


-* Xo 0 Xo -^*lU (23) 

i=1 


M d = -k Ro e Ro - fcn 0 en 0 + (Rq R 0d fl 0d ) A J 0 Rq R 0d fl 0d 


n 

+ J 0 RoR 0 po d - $r o 0 Ro - f' piRQ^Oxi, (24) 


for positive constants k XQ , k Xo , kR 0 , kn 0 G M. Here, the esti¬ 
mates of the unknown parameters 0 Xo , 0 Xi , 0r o are denoted 
by 6 Xq , Q x . , 0r o g W 10 . Adaptive control laws to update the 
estimates of disturbances are introduced later at Section IIII-I 
□ 


These are the ideal resultant force and moment to achieve 
the control objectives. One may try to choose the virtual 
control input pi by making the expressions in the right- 


hand sides of ( |2l| and (22), namely YU Pi and Y% PiRo Pu 
become identical to Fd and Md , respectively. But, this is not 


B. Design of Normal Components 

Substituting © into 0 and using the equation of 
motion for the i-link is given by 

A* = jqidi - f-qi{wi + ). (29) 

Here, the normal component of the control input uf is 
chosen such that qi qi d as t oc. Control systems for 
the unit-vectors on the two-sphere have been studied in [10], 
[11]. In this paper, we adopt the control system developed 
in terms of the angular velocity in [11], and we augment it 
with an adaptive control term to handle the disturbance A ^.. 

For the given desired direction of each link, its desired 
angular velocity is obtained from the kinematics equation as 

Ui d = Qi d x q id . 

Define the direction and the angular velocity tracking error 
vectors for the i-th link, namely e qi , e CJi G M 3 as 

e Qi = Qid x Qi, 
e Mi = w* + qpi d . 

For positive constants k q , k u G M, the normal component of 
the control input is chosen as 

uj- = mikqi{-k q e qi - k^e^ - ( q t ■ u> id )qi - qpd} 














- nrii^ai - ® X .0 X 


(30) 


Note that the expression of uf- is perpendicular to qi by 
definition. Substituting ( [30] ) into ( [29] ), and rearranging by 
the facts that the matrix — qf corresponds to the orthogonal 
projection to the plane normal to qi and qf = —q i9 we obtain 


uji = -k n e nA - k UJ e UJi - • o; id )^ - q*u d 


v q^qi 

1 

rrijl 




(31) 


where the estimation error is defined as Of. = 0 r 


0 X 


In short, the control force for the simplified dynamic 
model is given by 


Ui = u" + Ua 


(32) 


C. Design of Adaptive Law 

Next, we design the adaptive laws to construct the es¬ 
timates of unknown parameters. The following projection 
operator is introduced such that the estimated parameters stay 
in the bound of the true parameters given by ( p~8j ). 

V if ||0|| < B e 

Pr( 6 >, y) = { or ||0|| = B e and 0 T y < 0, 

(In e xn e ~ mW Sd T )y otherwise. 


Using this, the adaptive laws are defined as 

(33) 

Ox 0 = P r (^c 0 : y Xo ) 1 

(34) 

0r o = Pr {Or 0 ,Vr 0 ), 

(35) 

@xi — 

where y X0 ,yR 0 ,y Xi G R rae are defined as 

(36) 

y xo = + c x e Xo ), 

777/0 

(37) 

VRo = hRo®R 0 ( e no + c R e Ro)i 

Uxi = hxi*&x 0 \QiQi { (&XO + C x e xo) 

777/0 

(38) 

RoPi{en 0 + CReR 0 )} + <Ji( e ^ 

TTliii 

+ c q e Qi)] 5 


values of controller gains and controller parameters such that 
the following properties are satisfied. 

(i) The zero equilibrium of tracking errors 

(e Xo ,e Xo ,eR 0 ,e no ,e qi ,e Ui ) and the estimation errors 
{Oxq. 0r o , 6 Xi ) is stable in the sense of Lyapunov. 

(ii) The tracking errors asymptotically coverage to zero. 

(iii) The magnitude of the estimated parameters is less than 
Bq always, provided that the magnitude of their initial 
estimates is less than Be. 

Proof: Due to the page limit, the proof is relegated 
to [ 8 ]. ■ 

IV. Control System Design for Full Dynamic 
Model 

The control system designed at the previous section is 
based on a simplifying assumption that each quadrotor 
can generate a thrust along any arbitrary direction instan¬ 
taneously. However, the dynamics of quadrotor is under¬ 
actuated since the direction of the total thrust is always 
parallel to its third body-fixed axis, while the magnitude 
of the total thrust can be arbitrarily changed. This can be 
directly observed from the expression of the total thrust, 
Ui = —fiRie 3 , where fi is the total thrust magnitude, 
and R{e 3 corresponds to the direction of the third body- 
fixed axis. Whereas, the rotational attitude dynamics is fully 
actuated by the control moment M*. 

Based on these observations, the attitude of each quadrotor 
is controlled such that the third body-fixed axis becomes 
parallel to the direction of the ideal control force Ui designed 
in the previous section within a finite time. More explicitly, 
the desired attitude of each quadrotor is constructed as 
follows. The desired direction of the third body-fixed axis 
of the i- th quadrotor, namely b^ i G S 2 is given by 


63 , = - 


hi 


(40) 


(39) 

for positive constants c x ,CR,c q G M and adaptive gains 

O'xo 5 0jr q , h Xi G M. 

The first case of the projection map is the identity map, and 
the second case corresponds to the case that the estimated 
parameters are at the boundary of the region defined by 


This provides two-dimensional constraint on the three- 
dimensional desired attitude of each quadrotor, and there 
remains one degree of freedom. To resolve it, the desired 
direction of the first body-fixed axis b\. ( t ) G S 2 is introduced 
as a smooth function of time [7]. This corresponds to 
controlling the additional one dimensional yawing angle of 
each quadrotor. From these, the desired attitude of the i-th 
quadrotor is given by 


Bi c = 


(MX 




(18) and the unprojected direction y for the change of the 
estimates points outward. For such cases, y is projected onto 
the plane tangent to the boundary such that the estimated 
parameters remain on the region [ 12 ]. 

The resulting stability properties are summarized as fol¬ 
lows. 

Proposition 1: Consider the simplified dynamic model 
defined by M- For given tracking commands xo d , Ro d , a 
control input is designed as ([32])-([36]). Then, there exist the 


11*3*61* II 5 ' 

which is guaranteed to be an element of SO(3). The desired 
angular velocity is obtained from the attitude kinematics 
equation, Q ic = ( Rf c R ic ) w G M 3 . 

In the prior work described in [ 6 ], the attitude of each 
quadrotor is controlled such that the equilibrium Ri = 
Ri c becomes exponentially stable, and the stability of the 
combined full dynamic model is achieved via singular per¬ 
turbation theory [13]. However, we can not follow such ap¬ 
proach in this paper, as the presented adaptive control system 
guarantees only the asymptotical convergence of the tracking 
error variables due to the disturbances, thereby making is 










challenging to apply the singular perturbation theory. Here, 
we design the attitude controller of each quadrotor such that 
Ri becomes equal to Ri c within a finite time via finite-time 
stability theory [ 14 ], [ 15 ], [ 16 ]. 

Define the tracking error vectors e r. , en i G M 3 for the 
attitude and the angular velocity of the i-th quadrotor as 

e Ri = \( R I c R i ~ R I R i c ) V > e n, = - RfRi c Cl ic . 

2 

The time-derivative of e#. can be written as [7] 

= htr [ R J R i c ] I - Rf R ic )e Qi = E{R i ,R ic )en i . 

(41) 

For 0 < r < 1 , define S : M x M 3 -A M 3 as 

S( r iV) = [|j/i| r sgn(j/ 1 ), \y 2 \ r sgn(y 2 ), | 2 / 3 | r sgn(y 3 )] T , 

where y = [yi,y 2 il/ 3 ] T G M 3 , and sgn(-) denotes the sign 
function. For positive constants the terminal sliding 

surface G M 3 is designed as 

s% = eQ i + kRe r. + IrS(v , cr.). ( 42 ) 

We can show that when confined to the surface of Si = 0, 
the tracking errors become zero in a finite time. To reach 
the sliding surface, for positive constants k s ,l s , the control 
moment is designed as 



(a) 3D perspective 



(b) Top view 



Mi = —k s Si — l s S(r , Si) + x 

- ( k R Ji + IsrJidiagjUeR,, \ r ~ 1 ])E(R i ,R Ci )en i 

- JiihiRjR ic Q ic - RfRiAc)- ( 43 ) 

The thrust magnitude is chosen as the length of Ui, 
projected on to — it^e 3 , 


fi — 'U'i ‘ ^1^3') 


(44) 


which yields that the thrust of each quadrotor becomes equal 
to its desired value Ui when Ri = R ic . 

Stability of the corresponding controlled systems for the 
full dynamic model can be shown by using the fact that the 
full dynamic model becomes exactly same as the simplified 
dynamic model within a finite time. 

Proposition 2: Consider the full dynamic model defined 
by For given tracking commands xo d , Ro d and the 

desired direction of the first body-fixed axis &i., control 
inputs for quadrotors are designed as (43) and ( |44| ). Then, 
there exists controller parameters such that the tracking 
error variables (e Xo , e Xo , e Ro , e no , e qi , e Ui ) asymptotically 
converge to zero, and the estimation errors are uniformly 
bounded. 

Proof: See Appendix [A[ ■ 

This implies that the payload asymptotically follows any 
arbitrary desired trajectory both in translations and rotations 
in the presence of uncertainties. In contrast to the existing 
results in aerial transportation of a cable suspended load, 
it does not rely on any simplifying assumption that ig¬ 
nores the coupling between payload, cable, and quadrotors. 
Also, the presented global formulation on the nonlinear 


Fig. 2. Snapshots of controlled maneuver (redidesired trajectory, blueiactual 
trajectory). A short animation illustrating this maneuver is available at 
http ://youtu .be/nOWErfdzZLU 


configuration manifold avoids singularities and complexities 
that are inherently associated with local coordinates. As 
such, the presented control system is particularly useful for 
agile load transportation involving combined translational 
and rotational maneuvers of the payload in the presence of 
uncertainties. 


V. Numerical Example 

We consider a numerical example where three quadrotors 
(n = 3) transport a rectangular box along a figure-eight 
curve. More explicitly, the mass of the payload is mo = 
1.5 kg, and its length, width, and height are 1.0 m, 0.8 m, 
and 0.2 m, respectively. Mass properties of three quadrotors 
are identical, and they are given by 

rrii = 0.755 kg, ,J t = diag[0.0820, 0.0845, 0.1377] kgm 2 . 

The length of cable is U == 1 m, and they are attached to the 
following points of the payload. 

Pl = [0.5, 0, —0.1] T , 

pi = [-0.5, 0.4, -O.lf, p 3 = [-0.5, -0.4, -0.1] T m. 

In other words, the first link is attached to the center of the 
top, front edge, and the remaining two links are attached to 
the vertices of the top, rear edge (see Figure [I]). 

The desired trajectory of the payload is chosen as 

x o d (t ) = [1.2sin(0.27r£), 4.2cos( 0 .l 7 rt), — 0.5] T m. 
































(a) Position of payload (xo :blue, (b) Attitude tracking error of payload 
^o d :red) = §||#o - Ro d \\ 2 
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Fig. 3. Simulation results for tracking errors and control inputs, (for figures 
(c)-(f): i = l:blue, i = 2:green, i = 3:red) 


The desired attitude of the payload is chosen such that its 
first axis is tangent to the desired path, and the third axis is 
parallel to the direction of gravity, it is given by 


Ro d (t) = 


d 

11*0 J 


Initial conditions are chosen as 


g 3*0 d 
p3*o d || 



x o (0) = [1, 4.8, 0] T m, u o (0) = 0 3x i m/s, 

'/)(••) = e 3 , u>i( 0 ) = O 3 X 1 , Ri( 0 ) = / 3x3 , f 2 *( 0 ) = 0 3 x i. 

The uncertainties are specified as 

A* 0 = [1, 3, —2.5] t , Ar 0 = [-0.5, 0.1, —1.5] t , 
A Xi = [0.5, -0.2, 0.3] t , A Ri = [0.2, 0.3, -0.7] T . 

The corresponding simulation results are presented at 
Figures [2] and [3] Figure [2] illustrates the desired trajectory 
that is shaped like a figure-eight curve around two obstacles 
represented by cones, and the actual maneuver of the payload 
and quadrotors. Figure [3] shows tracking errors for the 
position and the attitude of the payload, tracking errors for 
the link directions and the attitude of quadrotors, as well as 


tension and control inputs. These illustrate excellent tracking 
performances of the proposed control system. 

Appendix 

a) Error Dynamics: From © and the dynamics 
of the position tracking error is given by 

n 

m 0 e Xo = m 0 {ge 3 - x 0d ) + A Xo + E^Aid + Af.). 

i= 1 

From ( [23] ) and ( [25] ), this can be rearranged as 

1 1 n 

= 9 e 3 — X 0d + + Y x + ——(A Xo + ^ 


m 0 


m 0 


i= 1 


1 


k Xo e x 0 kxo&xo + {*&x o 0x o + ^Xj^Xj) + Y x , 
m o ^ 

(45) 

where ( [13] ) is used and the estimation error is denoted by 
6 Xi = 0 X . — 0 X .. At the above equation, the last term Y x G M 3 
represents the error caused by the difference between q 7 and 
qi d , and it is given by 


1 n 

Y *~— E^f ~ V^icr 

m ° i =1 


We have /jLi d = q.i d q_l d H>i d from (28). Using this, the error 
term can be written in terms of e fl , as 


1 n 

Y * = ~ Qi d )Qi - QiJ 

m ° i=l 

1 " 

= E&iJ&e,.. 

0 • 


= 1 


Using ( [26] ), an upper bound of Y x can be obtained as 

1 n n 

\\y x \\ < — E ii^iiimi < E 7(11^11 + ii-v/iDik, i, 

iilo 


where 7 = 


i=1 


i=1 


. From (23) and (24), this can be 


m 0 y/\ m {VV T ]' 

further bounded by 


\\Y X \\ <Y^( k -o\\^ 0 \\+k ±0 \\e Xo \\) 

+ 7(^floll e «oll + feoll e n 0 ll) + -S}|l e gill) (46) 

where /3 = mo 7 , and the constant B is determined by the 
given desired trajectories of the payload, the assumption ( [17] ) 
on the bounds of <f> terms, and the adaptive law defined later 
that guarantee the boundedness of the estimated parameters 
0. Throughout the remaining parts of the proof, any bound 
that can be obtained from x 0d ,R 0d ,§Fj), or the adaptive law 
is denoted by B for simplicity. In short, the position tracking 
error dynamics of the payload can be written as ( [45] ), where 
the error term is bounded by <|46|». 

Similarly, we find the attitude tracking error dynamics for 
the payload as follows. Using ( [22] ), ( [24] ), and ( [27] ), the time- 
derivative of Joen 0 can be written as 

Jo^n 0 = (JqCq 0 + d) A en 0 - k Ro e Ro - k^e^ 

















































Define 


n 

+ *RoORo + ^ PiR 0 ^!r 0 ^» + (47) 

2=1 

where d = (2J 0 — tr[J 0 ] /)i?Q.^o d ^o d £ M 3 [9] that is 
bounded, and A# 0 G M 3 denotes the estimation error given 
by A r 0 = Ar 0 — A r 0 . The error term in the attitude 
dynamics of the payload, namely Yr G M 3 is given by 

n n 

Y R = Y1 ? iR 0 ( qtq i ~~ = - Y1 P iR 0 ( q LVi d )(li e <li- 

Similar with ( [46] ), an upper bound of Yr can be obtained as 

n 

||Yfl||<E^(^ol|exoll+fcxol|exoll) 

i=1 

+ Vi(kR 0 \\e Ro \\ + &n 0 ||en 0 ||) + -B}||e gi |], (48) 


where Si = mo ^ 
Next, from ( [3l| ), 


A m [?? T ] m o 


the time-derivative of the angular velocity 
error, projected on to the plane normal to q h is given as 


-qfe Ui =LO + (q • u d )q + q 2 cb d 

= —k q e qi — k UI e Ui — — — qi<$> x .0 Xi . (49) 

In summary, the error dynamics of the simplified dynamic 
model are given by (|47j) and ([49]). 

b) Stability Proof: Define an attitude configuration er¬ 
ror function for the payload as 

*Ro = |tr[/ - Ro d Ro] , 


which is positive-definite about Ro = Ro d , and ^r 0 = e# 0 • 
[7], [9]. We also introduce a configuration error function 
fy qi for each link that is positive-definite about qi = qi d as 




For positive constants e^ max , f>R 0 , f > qi , consider the following 
open domain containing the zero equilibrium of tracking 
error variables: 


77 - { (&X0 5 5 ^Ro "> ^0 5 5 5 5 ^-Ro 5 ) 

G (M 3 ) 4 x (M 3 x M 3 ) n x (M 3 ) 2 x M 3n | 

IK II < e x max j *R 0 < feo < 1, < K < !}• 

(50) 


In this domain, we have ||eR n || = \/\I/ rA2 — TK n ) < 
\ /fao(2-fa 0 ) = a 0 <l, and ||e 9i || = \/K( 2 “ K) < 
V^PJJ-RK) = oti < 1. It is assumed that f> Qi is 
sufficiently small such that naifd < 1. 

We can show that the configuration error functions are 
quadratic with respect to the error vectors in the sense that 


2 IKo l| 2 < ^Ro < 

glM 2 <K< 


1 


2 - Ko 
—-—I 

2 - K 1 


l e -Ro IP; 


where the upper bounds are satisfied only in the domain D. 


4^o— 2 ll^^o II 4" 2 ^oll e ^oll 4" c x e XQ • e XQ 

+ -en 0 • Jo^o + fcRo^o + • Joe^ 0 

n 1 

+ 2 ll e Kl| 2 + • e w ., 

2=1 Z 

where c x ,cr, c q are positive constants. This is composed of 
tracking error variables only, and we define another function 
for the estimation errors of the adaptive laws as 


V a = 


2 h r 


¥x 


2 h 


Ro 


HKii 2 +E^rii^n 


2=1 


The Lyapunov function for the complete simplified dynamic 
model is chosen as V = Vo + V a . 

Let = [|| e x 0 ||i ||ex 0 ll] T . = [|l e flol|»l|en 0 ||] T , % = 
[\\e q . II ? \\e COi ||] T G M 2 . The first part of the Lyapunov function 
Vo satisfies 


n 

zl 0 P X 0 Z X0 + Z% 0 P Ro Z Ro +y4^ Z 9. ^ V « 

2=1 

n 

< z T Xo P Xo z Xo + Z% 0 P Ro Z Ro + E4K-- 

2=1 


where the matrices P XQ , P Ro , P q ., P XQ , P Ro , P^. e M 2x2 
are given by 



kxo Cx 

-C X 1 J ’ 

2&#o_ -cr\ 
—cr\ X 



kxo Cx 
Y'x 1 

_ 2 Rrq_ 

2 -^h 0 

crX 


crX 

X 


R qi 


2 k„ 



1 

2 


2 kg 
2 ~^Qi 


c q 

1 


where A = A m [Jo] and A = Xm[Jo\- If the constants 
c x , cr 0 , c q are sufficiently small, all of the above matrices 
are positive-definite. As the second part of the Lyapunov 
function V a is already given as a quadratic form, it is 
straightforward to see that the complete Lyapunov function 
V is positive-definite and decrescent. 

The time-derivative of the Lyapunov function along the 
error dynamics ( [47] ), and ( [49] ) is given by 


17 — {k ±0 C x) ll^cco || Cx^Xo 11^0 II Cxfoxo^xo &Xo 

H - (Cx^xo H - &xo) ' Y x kflo II 4~ CR&R ’ Jo^-Qo 
— CRk Ro \\e Ro || 2 + c R e Ro • ((Joen 0 + d) A e^ 0 - fe 0 e^ 0 ) 
+ (e^ 0 + CrCr 0 ) • Yr 

n 

+ 'y ] ~{ku ~ Cg) || e UJi || — c q k q \\e q . || *— Cqk^e^ • e UJi 
2=1 

1 n 

+ —— (£®o + • ($*0^0 + ^ 

7710 2 = 1 

n 

+ (en 0 + c R e Ro ) ■ ($ Ro d Ro + y^pjPp 

i =1 



























“ ^2( e ^i+ C 1 e Qi) ' 

\i =1 


-^o-4 0 -E^'4, 

hR ° i =i 


mil, 

1 

l^Xn 


■**A 4 


h Oxo ' Oxo 

I^Xn 


Oxo ■ (yxo - Pr(^o> y*o)) + • (yflo - P^R 0 ,y Ro )) 

n 

+ 'y ] @xj ■ (y*i — Pr(0 ;Ei , y Xi ))- 




(O-O)-(y- Pt( 9, y)) = -—2(9-0). Wy 

Ini 

= j^(e T e-e T e)(e T y )< o. 

for each estimated parameter. The last inequality is due to 
(0 T 0 — 0 T 6) < 0 and (0 T y) > 0 obtained by (18) and d33]) . 

An upper bound of the remaining expressions of V at (51) 
can be obtained as follows. Since ||e# 0 || < 1, ||e^ 0 || < 
Il e n 0 || and ||d|| < B, 


V < -(k 


X 0 C x)\\ e x 0 


' k XQ ||c^ 0 1| C X k XQ C XQ ’ &x 


H - (cx^xo H - ^xo) ' Y x (kfi ,o 2cRA)||e^ 0 1| 
- CR&R 0 ||eR 0 || 2 + c R (kn 0 + ^)||e j R 0 ||||en 0 | 
+ {fn 0 + c R e Ro ) • Yr 




| C g^<7ll e Qill c qku e qi ' e Ui • 


i=l 


Cx^-Xq "b ^CCq) ‘ ^cll — 


Y.onPjcxk.. 

i =1 


N°iCo I 


c x^x 0 \\ e x 0 II ll^ico II + feo II^Xo II ) 


where z 7 ; = 
Wi G 


°XQ II 5 


ll^-Ro II> 


and the matrix 


(51) 


Wi = 


In the above equation, the expressions at the last four lines 
depending on the estimate error can be rearranged by using 
the adaptive laws, ([34])-([39]) as 


x3 is defined as 

^ \m\W Xi ] 

-h\\w xRt \ 

-hw xqi \\ 


-U W xRi\\ 

Xm[W Rz ] 

-UWnq.W 

where the sub-matrices are given by 
1 


AIM 

-|ll^ 9i || 

A m[W qi ] 


(56) 


W x . = 


W Ri = - 
n 


C-X^Xq (1 


noLif3) 

(1 + nai/3) 


C-x fox 


-(1 + noiifi) 


2 ^ -l i / outi j k x q (1 TIOLi /5) C, 

c R k Ro ( 1 - naidi) (/fa 0 + B + nc^cr*) 

- {kn 0 + B + (1 ~ ^cr*) - 2 crA 


From the definition of the projection map, the above expres¬ 
sions vanish for the first case of ([33]). For the second case, 


W qt = 


Cqkq 
Cq foou 
2~ 


Cq 

T~ 


km c a 


W xRi a 


W xqi = 


'J^x^Ro 3“ fiiC R k Xo 
lfk Ro + SiC R k Xo 
c x B 0 

e T _ -\- B 0 


lf c xkn 0 -b $ik x o 
jkn 0 + Sik Xo \ 7 
c r B 

&o <J ik Ro + B 


W xRi = 


Pkx 0 & I m ax 

If the constants c x , c R , c q that are independent of the con¬ 
trol input are sufficiently small, the matrices W Xi , W Ri , W qi 
are positive-definite. Also, if the error in the direction 
of the link is sufficiently small relative to the desired 
trajectory, we can choose the controller gains such that 
the matrix W{ is positive-definite, which follows that the 
zero equilibrium of tracking errors is stable in the sense 
of Lyapunov, and all of the tracking error variables Zi 
and the estimation error variables are uniformly bounded, 
i.e., g-xq 5 &xq 5 ^ Rq i £n 0 •> ^qi •> oi ? ^xq •> ^Rq •> ^Xi G ^oo’ These 


(52) 


From (46), an upper bound of the fourth term of the right- 
hand side is given by 


+ { c xB\\e x || + (Pk x0 e Xinax + B)\\e Xo ||}||e g< || 

+ ota(c x \\e xo \\ + ||e a:o ||)(fe jRo ||e J R 0 1| + ||e no ||). (53) 

Similarly, using ( |48| >, 

II (c R e Ro + en 0 ) ■ Y R \\ < 


also imply that e Xo , e Xo , e Ro , e no , e qi , e Ui € C 2 from ( |55| ), 
and that e Xo ,e Xo ,e Ro ,en 0 ,e qi ,e Ui e C^. According to 
Barbalat’s lemma [12], all of the tracking error variables 
e XQ , e Xo , e Ro , e ^ 0 , e qi , e w . and their time-derivatives asymp¬ 
totically converge to zero. 

A. Proof of Proposition [2] 

We first show that the attitude of the i-th quadrotor 
becomes exactly equal to its desired value within a finite 
time, i.e., Ri(t) = Ri c (t ) for any t > T for some T > 0. 
This is achieved by finite-time stability theory [14]. This 
proof is composed of two parts: (i) Si(t) = 0 for any t > T s 
for some T s < oo; (ii) when the state is confined to the 
surface defined by Si = 0 , we have e Ri (t) = e^ i (t) = 0 for 
any t > T R for some T R < oo. From now on, we drop the 
subscript i for simplicity, as the subsequent development is 
identical for all quadrotors. 

From [7], the error dynamics for eo is given by 


Y a i a i(c R k Ro \\e R J 2 +c R k n Je R J\\en 0 \\+kn 0 \\en 0 \\ 2 ) Je ^~ ^ x ^ + Af + A R + J{ClR T R c Vt c R T R c Cl c ). 

Substituting ( [43] ), 

§ 

Jen = —k s s - l s S(r, s) + A R - B s -t—t 


i= 1 


+ {c,R-B||e_R 0 || + (a 0 <Jik Ro + B)\\en 0 ||}||e 9i || 

+ a.A(c R \\e Ro \\ + \\en 0 \\)(k Xo \\e Xo \\ + k ±0 \\e X0 \\). (54) 

Substituting these into ([52]) and rearranging, 


v<Y~ z i w ^ 


~{k R J + l s rJ diag i [|e fl J r ])E(R,R c )e n . (57) 
Let a Lyapunov function be 
... 1 


(55) 


























From ( [42] ) and (J4TJ, its time-derivative is given by 
w = S • {Jen + (k R J + l s rJdmg j [\e Rj \ r - 1 })E(R,R c )e n }. 
Substituting ( |57) > and ( |43l >, and using (??), it reduces to 

>V = s • {-k s s - l a S(r, s) + A r - p||. B s } 

n 

<-fc s || S |i 2 -z 5 ^| Sj r +i +J B 5 || S ||-B 5 || S |j 

3 = 1 


where the last inequality is obtained from the fact that 


Fll“ < E"=i \ x i\ a for an y x = [xi,---, 
a < 2 [16, Lemma 2]. Therefore, 

>V < -eiW-e 2 W (r+1)/2 , 


and 0 < 


where 61 = iff/J wd e 2 = ^(A^ 7 [ 7 j) (r+1)/2 - This implies 
that s(t) = 0 for any t > T s , where the settling time T s 
satisfies 


T s < 


2 , ei W(0)i 1_r l/ 2 + e 2 

7^rj lU -«- 


according to [15, Remark 2]. 

Next, consider the second part of the proof when s = 
0. Let a configuration error function for the attitude of a 
quadrotor be 


J’ R = l -ix\l-R T c R]. 


Consider a domain give by Dr = {(R,Q) G SO(3) x 
R 3 \ ^/r < 'ipR < 2}. It has been shown that the following 
inequality is satisfied in the domain, 

\\\e R f<^R<D—\\e R \\\ (58) 

Therefore, it is positive-definite about e# = 0. The time- 
derivative of V R is given by 4 /r = 6r • e^. Therefore, when 
5 = 0, we have 


^r = -k R \\e R \\ 2 

< -&r||cr || 2 


Ir'E i W +1 

3 = 1 

^||efl|r +1 , 


Substituting ( [58] ), we obtain 

^R<-e^ R -e^ +1)/2 , 

where e 3 = 53 ^ and e 4 = ( 2 _^)(r + i )/2 • This implies that 
e R (t) = en(t) = 0 for any t > Tr , where the settling time 
Tr satisfies 

rr / 2 e 3 T'i?( 0 )( 1 - r )/ 2 +e 4 

Tr < —- -T In-. 

e 3 (l-r) e 4 

In summary, whenever t > T* = max{T s ,TR}, it is 
guaranteed that Ri(t) = Ri c (t ) for the i-th quadrotor. Next, 
we consider the reduced system , which corresponds to the 
dynamics of the payload and the rotational dynamics of the 


links when Ri(t) = Ri c (t). From (44) and (40), the control 
force of quadrotors when Ri = Ri c is given by 

~fi ■ Riez = (Ui ■ R Ci e 3 )R Ci e 3 = (■ -2G) - tt% = u, 


Therefore, the reduced system is given by the controlled 
dynamics of the simplified model. 

If the controller gains /c#, Ir, k s ,l s are selected large such 
that T* is sufficiently small, the solution stays inside of 
the domain D , where the stability results of Proposition 1 
hold, during 0 < t < T*. After t > T*, the controlled 
system corresponds to the controlled system of the simplified 
dynamic model, and from Proposition [I] the tracking errors 
asymptotically coverage to zero, and the estimation error are 
uniformly bounded. 
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